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Abstract. Models of the accretion disks of Young Stellar Objects show that they should not be ionized at a 
few AU from the star, and thus not subject to the MHD turbulence believed to cause accretion. This has been 
suggested to create a 'Dead Zone' where accretion remains unexplained. Here we show that the existence of the 
Dead Zone self-consistently creates a density profile favorable to the Rossby Wave Instability of Lovelace et al. 
(1999). This instability will create and sustain Rossby vortices in the disk which could lead to enhanced planet 
formation. 
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1. Introduction 

The existence of a Dead Zone in the accretion disk 
of Young Stellar Objects (YSO) was first proposed by 
|Gammie, 1996| he found that the inner region of the 
disk is ionized by collisions, and the outer region by 
cosmic rays. On the other hand there exists an in- 
termediate region, dubbed a Dead Zone, where the 
ionization remains so low that the gas is not cou- 
pled to magnetic fields. This creates a problem since 
accretion is widely believed to result from MHD in- 
stabilities, such as the Magneto-Rotational Instability 
( |Balbus fc Hawley, 1991| ). The Dead Zone should thus 
inhibit accretion in a significant portion of the disk. 
More recent stud ies (for example |D'Alessio et al., 1 998 
|Reyes-Ruiz, 2001| |Matsumura fc Pixdritz, 20050 have put 
constraints on the position and extent of the Dead Zone 
depending on the heating mechanisms and disk proper- 
ties, and related them to observations. Various estimates 
have been given for the extent of the Dead Zone, which 
is typically found to range between 1 and 5 AU from the 
star. 

In the present paper we discuss the consequences on 
the disk dynamics of the presence of such a zone, where 
the transport of matter and angular momentum (which 
we will represent by a turbulent viscosity) is significantly 
lower than elsewhere in the disk. In particular we show 
that an important role could be played by the Rossby 
Wave Instability (RWI), introduced by Lovelace et al. 



(1999) as a possible means of transporting angular mo- 
mentum in non-magnetized disks. This instability requires 
locally an extremum of a certain quantity, which we will 
discuss below. 

Here we show that the existence of a Dead 
Zone naturally generates such an extremum at its 
boundaries. The RWI naturally generates Rossby vor- 
tices. We relate this to the observation by sev- 
eral authors fBarge fc Sommeria, 1995| |Chavanis7"2 000 
|Johansen et al., 2004| ), that vortices can speed up planet 
formation, by allowing planetcsimals to accumulate in 
their center without following the viscous accretion. As 
shown in |Barge fc Sommeria, 1995] this might solve a 
problem encountered in present views of planet forma- 
tion, where planetesimals accrete to the central star on a 
viscous time scale, before they have had time to form a 
planet. We also show that the instability extends across 
the Dead Zone as spiral waves, which can maintain a sig- 
nificant accretion rate despite the absence of locally gen- 
erated turbulence. 



2. Rossby Waves Instability 

In 1999 Lovelace et al. (see als o Li et al. 2000, 2001), elab- 
orating on previous work by |Lovelace fc Hohlfeld, 197"8"1 
found a non-axisymmetric instability (Rossby Wave 
Instability, or RWI) occurring in unmagnetized disks, 
when there is a local extremum in the radial profile of 
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a quantity £. 
£ = T{r)S 2/T {r) 



Eft p 



(1) 



Where E is the surface density, the rotation fre- 
quency, k the epicyclic frequency and T is the 
adiabatic index. With T = 1, £ reduces to the 
inverse of the vortensity responsible for the corota- 
tion resonance of the Papaloizou-Pringle instability 
|Papaloizou fe Pringle, 1985) |Goldreich et al., 1986| 
Papaloiz ou*7TTin7^989y Numerical simulations by 
Ei^tHd^^OOOf confirmed that a bump or jump in the sur- 
face density E or the pressure p can cause the instability, 
and studied its consequences. Here we show how such a 
bump in the radial density profile can naturally result 
from the presence of a Dead Zone in the disk of a young 
stellar object (YSO). 



3. Numerical Model of the Dead Zone/ Instability 
Criterion 

3.1. Modelling the Dead Zone 

The Dead Zone is a limited region in the disk where the 
gas is not well coupled with the magnetic field, because 
the ionization is too low. If we make the assumption, as 
in |Gammie, 1996 ] that turbulent accretion is due to MHD 
instabilities such as the Magneto-Rotationnal Instability 
(MRI, Bal bus fc Hawley, 1991| ), we can model that region 
by a drop of the turbulent (MRI-generated) transport. 

We represent accretion in the disk by a turbulent vis- 
cosity (Sha kura fc Sunyaev, 19731 ) which drops in the re- 
gion of the dead zone, where MHD instabilities are in- 
active. We include this viscosity profile in the 2D non- 
linear hydrodynamic code developed by |Masset, 2 002 
|Masset fc Papaloizou, 2003] This is an Eulerian polar grid 
code with a staggered mesh and an artificial second or- 
der viscous pressure to stabilize the shocks (see also 
|Stone fc Norman, 1992"| ). The code uses the FARGO al- 
gorithm which significantly accelerates computations in a 
differentially rotating disk. The simulation are done us- 
ing standard boundary conditions, namely accretion is al- 
lowed through the inner boundary by using open bound- 
ary conditions and we assume the disk to be larger by 
allowing mass to flow in from the outer boundary. An ax- 
isymmetric, steady-state solution, which implies v E = 
Const., would give an artificially huge overdensity in the 
Dead Zone. Since the instability would cancel this condi- 
tion, we rather start from a simple power-law density pro- 
file that is allowed to evolve self-consistently during the 
simulations. Reaching a true steady state in these condi- 
tions would require prohibitive computation time and we 
rather concentrate in this paper on the development of the 
instability and a discussion of its consequences. 

The runs we present here are done with a resolution of 
N r = 150 and N$ = 450 grid points. We performed two 
types of simulations: in the first one, which allows only a 
limited resolution, the physical size of the grid spans radii 



between .25 AU and 10 AU, and the Dead Zone extends 
from 1 to 5 AU. In the second one, for higher resolution, 
and since we are mainly interested in what happens at the 
boundaries of the Dead Zone, we use a grid spanning from 
.25 to 4 AU, with a Dead Zone extending from .5 to 1 AU. 

In both cases the disk aspect ratio H/r = 0.04 where 
H . the vertical scale height of the disk, is uniform and 
constant. The sound speed of the gas is set from the 
disk aspect ratio. The initial density profile is E(r) = 
10 -5 (r/R m in) where R m in is the radius of the inner 
edge of the grid. Matter is allowed to flow in from the 
outer boundary, by accretion from the outer regions of 
the disk not included in the simulation. 

In order to model the Dead Zone we use an arctangent 
drop in the radial viscosity profile: 



v = v (e + atan(<5 r (r - r^f ))) if r < r, 



if r 



DZ 



v Q (e + atan(5 r (-r + r® n z ))) if r < r, 



DZ 
out 

> r > r 

DZ 



DZ 



where e represents a residual viscosity in the Dead Zone, 
and 5 r defines the steepness of the viscosity profile at the 
edges of the Dead Zone. All the runs presented here are 
done with e and S r equal to 10~ 5 and 50 respectively, 
giving the profiles shown on figure We did runs with e 
varying from 10 -2 to 10 -5 to confirm the behaviour with 
a less steep Dead- Zone. 







egular accretion 






regular accretion 


slow accretion zone | 





Fig. 1. Profile of the a- viscosity implemented to represent 
a Dead Zone between 1 and 5 AU with (e, S r ) = (10 -5 , 50). 



3.2. Evolution of the RWI Criterion 

With the prescribed viscosity profile, viscous accretion oc- 
curs in the outer region of the simulation but is far slower 
in the Dead Zone. As a result, gas accumulates near the 
outer edge of this zone, and naturally produces in less 
than a hundred orbits at that radius the extremum of £ 
required by the RWI (see fig|2]of the low resolution run) 
If we now look at the inner edge of the Dead Zone, the 
opposite happens, namely the matter in the Dead-Zone 
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is accreted slower than that in the inner region, therefore 
creating a dip inward from the Dead-Zone, also producing 
a favorable profile of £ 




Fig. 2. Left: plot of C a t = and once the vortices have 
formed, right: plot of the density in the disk at t= 100 
years, with the velocity field superimposed. The velocity 
field shows the vortex at 5 AU. There is also one at 1AU, 
not visible here. 



3.3. Formation of Vortices 

Once the threshold for the instability is reached we see 
Rossby vortices developping in the disk (there was still 
enough of the initial random perturbation to seed the in- 
stability). 



Figure shows a surface density plot of the disk at t = 
0, 100, 200, 300 years. We first see three vortices forming, 
but later only two of them surviving. These vortices are 



anticyclonic, and it 


is known (Barge & Sommcria, 1995 


Chavanis, 20001 


de la Fuente Marcos & Barge, 2001 


Johansen et al., 2004 


ll that this can cause an accumula- 



tion of dust and enhanced formation of planctesimals. 
Bar ge &: Sommeria, 1995| computed the amount of mass 
captured by a vortex depending on its distance to the 
star. Using our Dead Zone extending from 1 to 5 AU 
we see that the vortex at the inner edge of the Dead 
Zone will be able to capture 0.6 M ffi in 500 years (= 500 
orbits). From the same calculation the vortex at the outer 
edge of the Dead Zone will accumulate 16 M s in 6 10 3 
years (— 500 orbits). Those values, obtained in the case 
of a standard model nebula ( |Cuzzi et al., 19931 ) an d using 
dust simulations, are given here as an indication for the 
timescale/captured mass to expect from such vortices. 

3.4. Accretion Through the Dead Zone 

These simulations do not allow to follow the evolution of 
the vortices with enough precision. Therefore we focus our 
simulation box on either of the boundaries by putting it 
at the radius of r = 1 AU. This allows us to have a much 
better resolution at the location of the vortices. Figure 0] 
shows the evolution of the density profile in the region 
around the outer edge of the Dead Zone between t = 100 
years, when the RWI appears, and t — 400 years. This 




Fig. 3. Zoom of the first 2 inner AU of the simulation at 
t = 0,100,200,300 years, showing the density. One sees 
three vortices forming, later evolving to two vortices, near 
the outer edge of the Dead Zone. 



shows that accretion proceeds across the Dead Zone, de- 
spite the low viscosity. This can be attributed to the spi- 
ral waves generated by the vortices, and seen in figure 
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Fig. 4. This figure shows the density profile at t — 100 
(formation of the vortices) and 400 years : this shows that 
the spiral waves, seen in figure 2 to extend from the vor- 
tices, cause accretion across the Dead Zone. 
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J3J to extend across the Dead Zone. The generation of 
these waves corresponds to an exchange of energy and an- 
gular momentum between the vortices and the gas of the 
Dead Zone. We did not attempt to quantify this accretion, 
since this would require reaching a steady state which we 
can not obtain reliably with our present simulations, for 
a number of reasons: first, because of the 2D nature of 
the code we use; second, because on longer timescales the 
dynamics of the vortices will be influenced by the accu- 
mulated dust, and a two-phase code following separately 
both gas and dust would be required; finally, we also note 
that even a weak self-gravity would significantly enhance 
the coupling between the vortices and the spiral waves. 
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4. Conclusion 

Although the Dead Zone is not turbulent, we have shown 
that it naturally results in conditions supporting the de- 
velopment of the Rossby Wave Instability. The differential 
mass accretion creates a high density bump/ring at the 
boundaries of the Dead Zone and trigger the RWI.The vor- 
tices, in turn, survive long enough to be significant for dust 
accretion and enhanced planet formation. They also gen- 
erate spiral waves, resulting in significant accretion across 
the Dead Zone. 

Further studies will be needed to quantify these effects. 
This will require, in particular, a gas/dust /planetesimal 
code able to run on long time scales, and probably also 
the inclusion of self-gravity. 
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